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We experimentally demonstrate how coherent population trapping (CPT) for donor-bonnd elec¬ 
tron spins in GaAs results in autonomous feedback that prepares stabilized states for the spin 
polarization of nnclei aronnd the electrons. CPT was realized by excitation with two lasers to a 
bonnd-exciton state. Transmission studies of the spectral CPT feature on an ensemble of electrons 
directly reveal the statistical distribution of prepared nnclear spin states. Tuning the laser driving 
from blue to red detuned drives a transition from one to two stable states. Our results have im¬ 
portance for ongoing research on schemes for dynamic nuclear spin polarization, the central spin 
problem and control of spin coherence. 

PACS numbers: 03.65.Yz, 42.50.Ex, 76.30.Mi, 76.70.Hb, 78.47.jh 


Following the emergence of electron spins in quantum 
dots and solid state defects as candidates for spin qubits 
it has become a major goal to realize control over the nu¬ 
clear spins in such nanostructures. In many experimen¬ 
tal settings, interaction with disordered nuclear spins in 
the crystal environment is detrimental to the coherent 
evolution of carefully prepared electron spin states m- 
Preparation of nuclear spins in a state that has reduced 
spin fluctuations with respect to the thermal equilibrium 
state will help to overcome this problem [3]. Proposals 
to achieve this goal have been put forward for electron 
spin resonance (ESR) on one- or two-electron quantum 
dots 0111 , and for optical preparation techniques that 
either rely on a quantum measurement technique 00 
or a stochastic approach PE]. Experimental advances 
have been made with ESR and optical techniques on sin¬ 
gle quantum dots [T2UI8] and nitrogen-vacancy centers 
m, and on quantum dot ensembles [20l [2T] . 

Several of these works lEHni [13 [13 uni make use 
of the optical response of the electronic system near 
the coherent-population-trapping resonance (CPT, ex¬ 
plained below) because it is highly sensitive to perturba¬ 
tions from nuclear spins. Notably, these experiments so 
far have focussed on quantum dots where, due to the par¬ 
ticular anisotropic confinement, hyperfine coupling with 
a hole-spin in the excited state is reported to dominate 
m In recent work [31] we discussed how the interplay 
between electron-nuclear spin interaction and CPT influ¬ 
ences the stochastics of the nuclear spin bath for a class 
of systems where hyperfine interaction with the ground- 
state electron spin dominates. 

Here we report experiments on this latter class of sys¬ 
tems. We demonstrate an all-optical technique that sta¬ 
bilizes the nuclear spin bath around localized donor elec¬ 
trons in GaAs into a non-thermal state under conditions 


of two-laser optical pumping. We show that the nuclear 
spin system is directed either towards a single stable state 
or (probabilistically) towards one of two stable states, de¬ 
pending on laser detuning from the excited state. Our 
results show how feedback control arises from the inter¬ 
play between CPT and dynamic nuclear spin polariza¬ 
tion (DNP), and confirm that the electron-spin hyper¬ 
fine interaction dominates for our system (despite the 
strong similarity with the negatively charged quantum 
dot). Our results indicate that this interplay can be used 
to create stable states of nuclear polarization with re¬ 
duced fluctuations. 

We perform measurements on the nuclear spin dynam¬ 
ics in a lO-^m thick MBE-grown film of GaAs doped 
with Si donors at a concentration of ~3xl0^^ cm“^, 
which is well below the metal-insulator transition (at 
cm“^). The wafer is cleaved in 2-by-2 mm^ parts 
along the (110) crystal axes. The film is removed from 
a GaAs substrate by wet etching an AlAs buffer layer in 
HE. The film is then transferred to a sapphire substrate 
which allows us to do transmission measurements in a 
cryogenic microscope [23j . Measurements are performed 
at a temperature of T = 4.2 K and magnetic field of 
Bext = 5.9 T. The sample is mounted such that the mag¬ 
netic field direction is along the (110) axis. Light from 
tunable continuous-wave lasers (Coherent MBR-110) is 
delivered to the sample by a polarization-maintaining 
fiber and passes through the sample along the (100) axis. 
Transmitted light is collected in a multimode fiber and 
detected by an avalanche photodiode outside the cryo¬ 
stat. For getting reproducible data it was essential to 
stabilize laser powers within 1% and laser frequency drift 
within 10 MHz. 

The optical transitions that we address are from the 
donor-bound electron spin states (|t), I/)) to a level of the 
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FIG. 1. Energy levels and feedback control scheme, (a) Thick 
black lines are the spin states |1), |2) and optically excited 
state 1 3). Fs, 7 s and Fs ,73 are spin and excited state decay 
and dephasing rates, respectively. Two lasers at frequencies 
uj\, L 1 J 2 couple to the system with Rabi strengths Hi and 172 , 
excited state detuning A, and Overhauser shift 5 (see main 
text). The energy splittings cJis, a ;23 and uoz are fixed at 
the values for 5 = 0 (ft omitted for brevity), (b) Measured 
CPT signature in the n-GaAs sample (here for L 02 = UJ 23 and 
(5 = 0). (c) Left panel: two distinct control regimes for nuclear 
spin control. Middle panel: optically-induced electron spin 
polarization ((S'z) — (S^)) as a function of Overhauser shift 
5, with lasers fixed at cji = U113 + A and UJ2 = <^23 + A, 
displays two stable states of the nuclear spin bath for A < 0 
(red dots) and one stable state for A > 0 (blue dot). Right 
panel: expected steady state Overhauser shift distributions. 
Calculations with parameters 73 = 10, F^ = 10“"^, 7 s = 10“®, 
f7i = 172 = 0.5, A = ±1 normalized to r 3 = 1 | 22 |. 


bound trion {\D'^X)), that consists of two electrons and 
one hole bound at the silicon donor. These three states 
form a A-type energy level configuration, further defined 
in Fig.j^a). The magnetic field is applied perpendicular 
to the light propagation direction (Voigt geometry) such 
that the optical transitions have polarization selection 
rules discriminating between horizontally (<7+,-, coupling 
to It) —|T*°V)) and vertically (tt, coupling to |t) —|Z7°V)) 
polarized light. 

We first demonstrate CPT for our system. CPT is 
a narrow resonance in two-laser driving as in Fig. [^a) 
[21] where the system gets trapped in a dark state (for 
ideal spin coherence I'F) oc O 2 |t) — |i))- In trans¬ 

mission this appears as a narrow window of increased 
transparency within the broader absorption dip when one 
laser is scanning while the other is fixed (Fig. [^b)). Its 
position signals two-photon resonance, and occurs where 


uji — UJ 2 equals the the electron spin splitting. The line- 
shape of the CPT resonance can reveal information about 
the electron spin states, which can be obtained by fitting 
the curve in Fig. [^b) to the Lindblad equation for the 
A-system [2S]. Our n-GaAs samples yield an inhomo¬ 
geneous dephasing time T| « 3 ns [21]. However, the 
homogeneous dephasing time T 2 has been estimated to 
be at least 7 /rs [26] with a spin-echo technique. The dis¬ 
crepancy between T 2 and T 2 is largely due to dephasing 
caused by ~ 10 ® disordered nuclear spins per electron. 

Due to the Fermi contact hyperfine interaction, a non¬ 
zero nuclear spin polarization exerts an effective magnetic 
(Overhauser) field on the electron spin and causes a 
shift of the electron spin levels, denoted by S in Fig. [^a). 
The value of (5 = p^max is proportional to the nuclear 
spin polarization p G [— 1 , 1 ], where (5max is the maximum 
shift set by the hyperfine interaction strength. For the 
donor electron in GaAs 6max = 24.5 GHz (obtained from 
the maximum Overhauser field [27j via S = 
with g-factor g = —0.41 [21] )• The thermal equilibrium 
properties of the nuclear spin bath are well approximated 
by considering N non-interacting spins / with gyromag- 
netic ratio 7 . Then p and its variance are in the 
high temperature limit hjBext/kBT <C 1 (our experi¬ 
mental conditions) p = Ti'^Bext{I + 1)/3A:bT « 0 and 

al = (/+l)/3/V-p2 g]. 

Because nuclear spin dynamics is slow as compared to 
the electron’s, light interacting with the system sees a 
snapshot of the Overhauser shift taken from a distribu¬ 
tion P{5). A measurement on an ensemble of these sys¬ 
tems should account for averaging over P{5). The CPT 
lineshape of Fig. [^b) arises from the transmittance, with 
a susceptibility that is averaged over P{5), 

/ /7 \ 

T{u:i) = exp S)dSj , (1) 

where d is the thickness of the medium, p the density of 
donors, c the speed of light, i = I, 2 labels the laser fields. 
Here Xi is the susceptibility for the laser field for a fixed 
(5. It can be calculated from the Lindblad equation and 
depends on other system parameters implicitly |22j . At 
thermal equilibrium P{5) is a Gaussian centered at zero 
with variance tr^ = ^max^p. For 1 = and N = 10® 
it has a width (FWHM) of 2-^2 log{2)as = 136 MHz, 
which roughly corresponds to the width of the measured 
CPT. 

However, P{6) can undergo changes when the elec¬ 
tron spin is brought out of thermal equilibrium by op¬ 
tical orientation. An optically-induced electron spin po¬ 
larization will in turn induce nuclear spin polarization 
via a hyperfine-mediated cross-relaxation process known 
as DNP. In Ref. [21] it was described how the interplay 
between the laser-induced electron spin polarization near 
CPT resonance and DNP can change the shape of P{S) 
by autonomous feedback control, leading to the forma¬ 
tion of stable states for the nuclear spin polarization 
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and offering the potential of reducing the variance cr|. 
The essence of this method is pictured schematically in 
Fig. [^c). It shows two distinct control regimes (color 
coded, red and blue) where both lasers are either red 
(A < 0) or blue (A > 0) detuned from the excited 
state. The change in laser coupling strength with S is 
asymmetric when A 7 ^ 0 (one laser approaches reso¬ 
nance while the other moves away from it). For a single 
system with a particular Overhauser shift this causes a 
sharp change in the optically-induced electron spin po¬ 
larization {Sz) — (Sz) (where the overbar implies that 
the expectation value is taken at thermal equilibrium), 
shown in the middle panels as a function of 6 (the Over¬ 
hauser shift is here normalized to Fa). The blue and 
red dots indicate stable points, where {Sz) = {Sz) and 
d/dS{{Sz) — {Sz)) < 0. We thus expect P{d) to evolve 
from the initial Gaussian to either a distribution with 
two maxima, or to a distribution with one maximum. 
Such steady-state distributions are non-thermal and can 
thus have reduced fluctuations if the system’s feedback 
response (slope of {Sz) — {Sz) near the stable point) is 
strong enough j^ . 

We investigate this interplay between CPT and DNP 
for the donor-bound electrons in GaAs by monitoring the 
changes in the CPT lineshape induced by two-laser op¬ 
tical pumping, with both lasers at equal intensity near 
two-photon resonance. Figures [^a,b) show the CPT 
lineshape before (gray lines) and after 10 min of optical 
pumping with blue- and red-detuned lasers. While scan¬ 
ning over the ensemble CPT peak, the probe laser meets 
exact two-photon resonances (near-ideal CPT peaks) of 
individual electrons for a range of d-values. The sus¬ 
ceptibility is thus proportional to the number of elec¬ 
tron spins experiencing a particular Overhauser shift S, 
hence reflecting the underlying nuclear spin distribution. 
The nuclear spin distribution stabilizes as predicted in 
both cases, observed as a non-shifted single CPT peak in 
Fig. da) and a non-shifted split CPT peak in Fig. [^b) 
(the splitting directly reflects the doubly peaked P{d) of 
Fig. 0 c)). This is in clear contrast with a CPT peak 
recorded after 10 min of single-laser optical pumping (in¬ 
set Fig. [^a)), which shifts the CPT peak by ~400 MHz 
since DNP gives here a net nuclear spin polarization. 

The lineshape in the main panel of Fig. [^a) remains 
similar, while a narrower and higher CPT peak is ex¬ 
pected if the width of the stabilized P{S) would indeed 
be reduced. In Ref. [52] it was pointed out that for an 
open system the narrowing by the feedback mechanism is 
in competition with nuclear spin diffusion. For donors in 
GaAs this plays a stronger role than for quantum dots, 
where a material barrier surrounding the dot suppresses 
this spin diffusion. Not observing a narrowing of the CPT 
peak is also due to non-uniform laser intensities for the 
electron ensemble (further discussed below). 

Figure [^c) shows the transition from red- to blue- 
detuned two-laser pumping, for a range of detunings A. 



Detiining A (GHz) Probe laser scan (MHz) 


FIG. 2. CPT signatures of DNP feedback, (a) The CPT 
peak in the transmission signal as a function of probe laser 
frequency, before (gray) and after (blue) a DNP pumping pe¬ 
riod with two lasers fixed on two-photon resonance (Fig. 1(c)) 
and detuning A = -1-4 GHz. The inset shows how DNP pump¬ 
ing by one laser (on the |t) — \D^X) transition) causes a shift 
of the CPT resonance. Two-laser pumping stabilizes the nu¬ 
clear spin polarization at its thermal equilibrium value (here 
without observing a significant narrowing), (b) Results ob¬ 
tained as for panel (a), but with A = —6 GHz. In this case 
the CPT peak after DNP pumping (red) shows a splitting, 
(c) CPT traces taken after DNP pumping, for various values 
of A. (d) Values of the peak splitting, obtained from traces 
as in panel (c). All data was taken with both laser intensities 
stabilized at values of about 3 Wcm“^ (for DNP pumping and 
CPT probing). Black line: simulation with parameters as in 
Ref. |22| . except Fd/F/j = 4000 and 73 = 20 GHz. 


Splittings in these CPT peaks are analyzed in Fig.j^d), 
obtained by fitting two Gaussians to each CPT peak. 
Where the fit does not improve with respect to a single- 
Gaussian fit we take the splitting to be zero. The data 
reproduces the essential features of the model [ 22 ] (black 
line), showing a discontinuous transition and a maximum 
splitting when the pump lasers are tuned to slope of the 
transition line at A « —5 GHz, where the response to 
a shift of 5 is largest. We analyzed that this transition 
is a unique feature that confirms the dominance of the 
electron spin for the relevant DNP mechanism [22] • For 
A > 0 there is no good match, but the fitting also yields 
larger error bars. We attribute this to inhomogeneous 
broadening in the optical transitions (effective spread in 
detunings A) which prevents all systems from making the 
transition simultaneously. 

We now focus on the control regime A < 0 to examine 
the dependence of the stabilization on the control pa- 
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Rabi frequency (GHz) 




FIG. 3. (a) CPT traces after DNP pumping with two lasers at 
A « —3.5 GHz. The trace labeled cui = 0 MHz is taken after 
pumping on exact two-photon resonance. The blue (red) trace 
is measured after pumping with uji = )31 MHz detuned 

from exact two-photon resonance (see Fig. 1(c)). (b) The 

peak splitting in GPT traces after pnmping with A « —4 GHz 
and exact two-photon resonance, as a function of the inten¬ 
sity of the two lasers (keeping the intensity ratio fixed near 
1). The gray background shows the range where the GPT 
peak shape was analyzed as a single peak. Significant double¬ 
peak character was observed for the total laser intensity above 
~3 Wcm“^. GPT traces were all taken with both lasers in¬ 
tensities at ~3 Wcm“^. Black line: simulation with same pa¬ 
rameters as in Fig. the top axis shows the Rabi frequency 
corresponding to the simulation [22j . 


rameters during the optical pumping phase. Figure]^ a) 
shows the importance of carefully tuning the relative fre¬ 
quencies for getting a balanced distribution. A detuning 
as small as 31 MHz for one of the lasers gives a signifi¬ 
cant shift within P{S) to either one of the stable states. 
Figure|^b) shows values for the splitting as a function of 
the laser powers (varied simultaneously). The splitting 
shows a discontinuous onset and subsequent increase due 
to power broadening of the CPT peak. The data quali¬ 
tatively matches the prediction ([H], black line) but the 
slope is lower than the simulation. We attribute this to 
standing wave patterns in the GaAs layer (which acts as 
a weak cavity). The patterns for the two lasers do not 
fully overlap since they differ in frequency. This prohibits 
addressing the entire ensemble with equal laser intensi¬ 
ties, and gives for the ensemble an averaged, less effective 
feedback mechanism. This also provides a limitation for 
the amount of CPT-peak narrowing in the blue-detuned 
case. The narrowing effect relies on carefully balanced 
laser intensities, and this is compromised due to the in¬ 
tensity variation inside the sample. Studying the achiev¬ 
able narrowing of P{6) requires an experiment with uni¬ 
form intensities for the ensemble. 

Figure presents time evolution of the effects. Fig¬ 
ure Qa) shows build-up of the splitting, obtained by tak¬ 
ing CPT traces during the optical pumping phase every 
30 s (each trace is collected within 1 s). The splitting 
stabilizes after approximately 4 min. Figure [^b) shows 
decay of the splitting. It consists of traces collected af¬ 



Probe laser scan (MHz) 


FIG. 4. Time evolution of build-up (a) and decay (b,c) of 
stabilized nuclear spin ensembles, measured during and after 
the DNP pumping period (subsequent traces top to bottom, 
as labeled). Panel (b) gives a reference for the GPT peak 
before pumping. The data in (a) is obtained from CPT scans 
of 1 s in between periods of 30 s DNP pumping with two lasers 
fixed at two-photon resonance (A = —4 GHz). Panel (b) 
presents CPT scans of 1 s, taken after a fully dark period of 
1 s (top trace) to 10 s (bottom trace) after DNP pumping 
(A = —4 GHz). The results in (c) are from continuously 
taking CPT scans of 1 s (only four traces shown), after a DNP 
pumping period at A = —2 GHz. All data was taken with 
both laser intensities stabilized at values of about 3 Wcm“^ 
(for DNP pumping and CPT probing). 


ter the optical pumping phase. After 10 min of optical 
pumping (repeated before each trace) the system is kept 
in the dark for a time ranging from 1 to 10 s. The split¬ 
ting fades away in seconds, consistent with the relaxation 
of the lattice nuclear spins by spin diffusion away from 
the electron [5S]. However, when CPT scans are taken 
continuously after the optical pumping phase the split¬ 
ting decays much slower and persists up to at least 20 min 
(Fig-i c), we verified that taking such scans without the 
preceding pumping phase does not induce a splitting). 
We attribute this to a suppression of the spin diffusion 
while the system is illuminated: under optical excitation 
(during CPT scans) the electron spin is most of the time 
significantly polarized and this suppresses nuclear spin 
diffusion because it creates an inhomogeneous Knight 
field for the surrounding nuclear spins [29l |30] . This ef¬ 
fect could be used to improve the strength of the feedback 
control and the amount of narrowing: if the temperature 
of the experiment would be lower or the magnetic field 
stronger (increased (Sz)) the thermal-equilibrium elec¬ 
tron spin polarization can suppress nuclear spin diffusion. 

Our results open the possibility to use the interplay 
between CPT and DNP to operate a mesoscopic spin 
system as a feedback loop that converges towards a well 
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defined steady state, determined by laser power and de¬ 
tuning, with the possibility of reduced nuclear spin fluc¬ 
tuations and less electron spin dephasing. The mecha¬ 
nism is generally applicable to localized spins where DNP 
is dominated by electron-nuclear spin hyperfine coupling 
and can also be used for other paramagnetic defects, as 
ensembles or single systems. A notable example is the flu¬ 
orine donor in ZnSe [HI |32] , a II-VI material with dilute 
nuclear spins (in GaAs all atoms have non-zero nuclear 
spin). Nuclear spin diffusion, mediated by dipole-dipole 
interaction (inversely proportional to distance between 
nuclear spins to the power 6), will here be much less a 
limitation for narrowing. 
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